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Particle acceleration close to the supermassive black hole horizon: the case of M87 
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The radio galaxy M87 has recently been found to be a rapidly variable TeV emitting source. We analyze the implica- 
tions of the observed TeV characteristics and show that it proves challenging to account for them within conventional 
acceleration and emission models. We discuss a new pulsar-type scenario for the origin of variable, very high energy 
(VHE) emission close to the central supermassive black hole and show that magneto-centrifugally accelerated electrons 
could efficiently Compton upscatter sub-mm ADAF disk photons to the TeV regime, leading to VHE characteristics 
close to the observed ones. This suggests, conversely, that VHE observations of highly under-luminous AGNs could 
provide an important diagnostic tool for probing the conditions prevalent in the inner accretion disk of these sources. 
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1. Introduction 

The advance in modern imaging atmospheric 
Cherenkov technologies has opened up a new win- 
dow on the non-thermal Universe, providing a unique 
diagnostic tool to study the physical conditions in 
violent regions of galactic and extragalactic sources. 
Until recently, blazars have been the only extragalac- 
tic sources known to emit TeV (10^^ eV) radiation. 
Although comprising only a small fraction of radio- 
loud active galactic nuclei (AGNs), the favourable 
orientation of their jets at small viewing angles i 
to the line of sight {i < 10°), combined with the 
bulk relativistic motion of the plasma in their jets, 
makes them a privileged source class on the TeV sky. 
The H.E.S.S. detection of variable TeV emission from 
the non-blazar (i > 19°) radio galaxy M87 (redshift 
z — 0.0043, distance d ~ 16 Mpc) during a campaign 
in 2003-2006 [1], may thus come surprising. Although 
some detectable TeV emission has been anticipated 
based on simple homogeneous SSC model extrapo- 
lations [2], the observed rapid 7-ray flux variations 
and the hard TeV spectrum comes quite unexpected. 
The giant elliptical galaxy M87 is known to host one 
of the most massive black holes Mbh — 3 x 10^ Mq 
in the Universe, and a prominent, non-aligned jet 
detectable from radio to X-ray wavelengths. As we 
show below, the VHE results may point to a new 
acceleration process occurring in the vicinity of the 
central supermassive black hole (BH). 



2. Observed VHE characteristics in 
M87 and implications 

The H.E.S.S. observations of M87 in 2005 [1] show 
a gamma-ray spectrum that reaches beyond 10 TeV 
and is consistent with a relatively hard power-law 
(with spectral index a 1.2, where Si, oc 
The detected VHE output is quite moderate, with an 
isotropic TeV luminosity of ^ 3 x 10''" erg/s. During 
the 2005 high state, significant variability (flux dou- 
bling) on time scales of A^obs ~ 2 days was found, 
the fastest variability observed in any waveband from 
M87 so far. 

Several scenarios for the origin of VHE emission from 
M87 have been discussed in the literature, ranging 
from the vicinity of the central supermassive BH [3, 
4] and the innermost part of the jet [5-7] to the su- 
perluminal feature HST-1 at a de-projected distance 
of 100 pc from the central BH [8, 9]. The latter 
has some interesting features including the fact that 
the TeV flux maximum in 2005 seems to coincide 
with a Chandra X-ray flux maximum from HST-1. 
On the other hand, the same observations [1, 10] 
also indicate that the TeV and HST-1 X-ray flux 
evolutions from 2003 to 2004 are not correlated, yet 
the TeV spectral shape does not change (i.e., same 
power index within errors) from 2004 to 2005, so 
that an HST-1 origin of the TeV emission in 2005 
seems to require some sort of cosmic conspiracy. 
Note further, that the observed spectrum of HST- 
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1 appears synchrotron cooling dominated above the 

break at ~ lO""^^ Hz, where the spectrum steepens 
to a ~ 1.0 {Si, oc [10], corresponding to an av- 

erage radiating particle distribution 71(7') oc 7'"^. In 
order to Compton upscatter lower energy photons to 
the observed TeV regime, the presence of energetic 
electrons with Lorentz factors up to 7^ ~ lO'^/S is re- 
quired within the source region, where 6 is the bulk 
Lorentz factor of the source. If scattering in a sin- 
gle SSC-type scenario occurs in the Thomson limit, 
the target photon field, and thus the observed TeV 
flux, is expected to vary on time scales larger than 
implied by a source limit inferred from the cooling 
break, i.e., the total source size would be at most 
comparable to i? ~ c Atobs ^- If) on the other hand, 
the relevant scattering occurs in the Klein-Nishima 
regime, the resultant TeV spectrum would be very 
steep, contrary to the observed one. Further studies 
are required to assess whether current HST-1 models 
are flexible enough to overcome these issues. 
A similar consideration seems to hold for leptonic 
scenarios that assume Compton scattering in a decel- 
erating relativistic inner jet flow [5]. Radio features 
in the inner jet, for example, seem to be at most 
mildly relativistic [11, 12], yet significantly relativis- 
tic on larger scale [8], the jet thus more resembling 
an accelerating than a decelerating flow, although 
matters may be more complicated if an internal ve- 
locity structure (e.g., spine plus sheath) is present. 
Note also, that the predicted spectral shapes for the 
TeV regime in simple model applications are much 
steeper than observed, although again, an internal 
jet structure may help by introducing additional ac- 
celeration and radiation effects [7, 13]. 
Here we consider an alternative scenario assuming 
particle acceleration and TeV 7-ray production to 
take place close to the event horizon of the central 
supermassive black hole. We show that such a model 
can successfully reproduce the observed VHE char- 
acteristics and introduce an important link between 
accretion disk physics and jet formation theory. 

3. Particle acceleration close to the 
centred black hole 

Today, magneto-hydrodynamical models are widely 
considered to represent the most promising class for 
the formation and coUimation of relativistic astro- 
physical jets. According to this picture, magnetic flux 
dragged inward and amplified by dynamo actions in 



the inner accretion disk can build up a rigidly rotat- 
ing, dipolar magnetosphere. Along open flux surfaces 
the bulk of the plasma is centrifugally accelerated to 
relativistic speeds {Tf, ~ 10) and collimated outside 
the Hght cylinder tl — (5 — 10) r^, e.g., sec [14, 15]. 
It has been realized for quite a while [16-18], that 
such a MHD field structure could also allow for ef- 
ficient centrifugal acceleration of test particles, i.e., 
a test particle co-rotating with the field line (bead- 
on-wire motion) will experience the centrifugal force 
and gain rotational energy while moving outwards 
[19, 20]. The radial motion of such a particle is most 
conveniently analyzed in the framework of Hamil- 
tonian dynamics [4, 17]. Consider an idealized two- 
dimensional model topology: since the Lagrangian L 
for a particle with rest mass mo on a relativistically, 
rigidly rotating wire (angular velocity Q = c/ri^ = 
constant) is not expHcitly time-dependent, i.e., 

L = -moc^^l - r'^/rl - jc? , (1) 

the associated Hamiltonian 

H = jnio c^(l - r'^/rl) = constant , (2) 

with 7 = (1— r^/r£— r^/c^)""'^/^ the Lorentz factor, is 
a constant of motion. Thus, as a particle approaches 
tl , the term in brackets of Eq. (2) gets smaller and 
has to be compensated by an increase in 7, which im- 
plies that the Lorentz factor increases dramatically 
for a particles approaching the light cylinder. In real- 
ity, unlimited growth will be prohibited by radiative 
energy losses (e.g., inverse Compton in the ambient 
photon field), the breakdown of the bead-on-the-wire 
(BW) approximation or the bending of the field line 
with increasing inertia [16, 17]. In particular, validity 
of the BW approximation requires, that the charac- 
teristic acceleration time scale, which can be easily 
derived from Eq. (2), i.e., 

is always larger than the inverse of the relativistic 
gyro- frequency co = eB/jnioC. Here m is deter- 
mined by the initial conditions, i.e. m = 1/(7q [1 — 
Tg/r^]^). This constrains achievable Lorentz factors 
to 

which is formally equivalent to the requirement that 
the Coriolis force must not exceed the Lorentz force 
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[17]. For parameter relevant to M87 (i.e., -B(rL) ~ 10 

G, TL - 5 X IQi'^ cm), one finds -/^^^ - 5 x 10^ 
for electrons. In general, centrifugal acceleration is 
less favourable for protons in the sense that 7max 

2/3 

ttiq ' indicating that accelerated protons will not be 
able to interact efficiently with the ambient photon 
field. 

4. Inverse Compton upscattering of 
Comptonized disk photons 

In realistic astrophysical environments, radiative en- 
ergy losses will always compete with acceleration 
and introduce additional constraints on the maxi- 
mum achievable particle energies. Given its very low 
bolometric luminosity output, mass accretion in M87 
is likely to occur in a two-temperature, advcction- 
dominated (ADAF) mode [14, 21, 22]. For the esti- 
mated Bondi accretion rate of to ~ 1.6 x 10~^TOEdd 
[21], the characteristic ADAF spectrum close to Tl 
peaks at aroimd i^'^ ^ 10^^ Hz with an associated 
(radio) luminosity Lr ~ 10^° erg/s. Inverse Comp- 
ton interactions with the ambient photon field are 
thus not expected to lead to a much stronger con- 
straint on the electron Lorentz factors compared with 
7max: suggesting that electron Lorentz factors up to 
7 ~ (10^ — 10^) could well be achieved via the pro- 
posed acceleration mechanism [4]. Hence, Compton 
(Thomson) upscattering of sub-mm {ly' < ~ 10^'^ 
Hz) accretion disk photons can easily result in VHE 
photons with energies extending up to ^ 10 TeV and 
beyond. Note that this is different from the case of 
luminous AGN sources, where severe inverse Comp- 
ton losses tend to limit achievable electron Lorentz 
factors to be below one thousand [16, 17]. Yet, even 
for the latter sources, centrifugal acceleration could 
operate as a promising pre-acceleration mechanism, 
possibly along with electrostatic wave surfing [23], 
and thus provide the energetic seed particles required 
for efficient Fermi acceleration on larger scales. 
For the inferred Bondi accretion rate rhs in M87, 
comptonization of cyclosynchrotron soft photons in 
an ADAF is expected to add a power law-like tail 
to the disk spectrum above u'^ with spectral index 
close to ac ~ 1-2 [4, 24]. Suppose that energetic 
electrons, after being released from the acceleration 
process with high Lorentz factors, encounter such a 
comptonized disk photon field. In general, the shape 
of the emergent (singly, Thomson) scattered inverse 
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Fig. 1. Differential energy spectrum of M87 for the H.E.S.S. 
2005 data. Error bars given arc for statistical errors only. 
Shown is a characteristic TcV spectrum in the range 47|j^( < 
u < 47I1/2, resulting from (Thomson) Compton upscatter- 
ing of a seed photon spectrum with spectral index ac = 1.2 
between i^J < i^' < 1^2, with i^(/i^2 ~ 0.015, p = 2 and 
72 ~ 2 X 10'''. 

Compton spectrum iic(z^) will depend on both the 
seed photon spectrum F{u') and the differential elec- 
tron distribution n{'y), e.g., 

jiciiy) OC r dv'F{v') i/'(p-i)/2 

/•minjxi ,1} 

X / dxx^^-^y^hc{x), (5) 

J min{x2 ,1} 

assuming a power law electron distribution 71(7) oc 
7~P in the interval 71 < 7 < 72 ~ 5 x 10^, with 
Xi := v/{4'yfiy') and ficix) = 2xlnx -|- x + 1 - 2x'^. 
Far away from the endpoints, i.e., for v <^ 47|i/^ ~ 5 
TeV, we obtain the common jic('^) oc i/~(p~i)/2 
power law evolution, independent of the incident 
photon spectrum. Yet, closer to the endpoints, i.e., 
for 4:^2^ < < 47|i/2, we become again sensitive 
to the seed photon distribution so that the resultant 
inverse Compton VHE spectrum will be power law- 
like with spectral index a ~ ac. This is illustrated in 
Fig. 1 for p = 2. 

5. Discussion 

Taking the detected TeV emission in M87 to arise via 
inverse Compton scattering by centrifugally acceler- 
ated electrons close to the central black hole, requires 
to examine whether, amongst others, the model can 
self-consistently account for the observed variability 
and luminosity output: 
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(1) The observed fliix doubling time scale is of or- 
der ^ 2 days. Centrifugal acceleration, on the other 
hand, occurs on a characteristic time scale tl/c. The 
observed variability thus implies a light cylinder ra- 
dius of order tl ~ 6 , which is indeed consistent 
with expectations from relativistic MHD jet forma- 
tion models [14, 15]. 

(2) The detected TeV luminosity is of order ixeV — 

3 X 10^° erg/s (if assumed isotropic). It can be shown 
[4] that the test particle number density He, required 
to account for ixeVi does not violate the presumed 
quasi force-free MHD condition in the sense that the 
corresponding kinetic energy density ~ ne^2^ec'^ is 
still much smaller than the allowed particle energy 
density of the background plasma. 

(3) Is it possible that TeV photons, produced via 
inverse Compton upscattering, can escape from 
the vicinity of the central black hole? In prin- 
ciple, photons of energy E [TeV] will interact 
most efficiently with target photons in the infrared 
regime, i.e., with energy e ~ (ITeV/E) eV. If 
the estimated, observed infrared luminosity Lir. ~ 
10**^ erg/s [25] is produced on characteristic scale 
iJiR, the optical depth for 77 absorption becomes 
r^^ ~ 0.2(LiR/104ierg/s)(ri/RiR)(E/lTeV), indi- 
cating that TeV photons can escape even if a non- 
negligible percentage of the observed infrared flux 
would originate in the innermost region of the ac- 
cretion flow ri ^ 60 Tg. In particular, assuming all 
of the observed flux in M87 to arise in an ADAF 
(which is perhaps over-restrictive given its jet) con- 
straints transparent radii of 10 TeV photons to scales 
between 5 and 13 Schwarzschild radii [26], consistent 
with the scenario presented here. 

6. Conclusion 

VHE gamma-ray observations of low-luminous, non- 
blazar AGN jet sources like M87 could provide an 
ideal test laboratory for particle acceleration close to 
the supermassive black hole horizon and serve as a 
fundamental diagnostic tool for the conditions preva- 
lent in the innermost part of the accretion flow. In 
blazars, these effects are likely to be swamped by rel- 
ativistically beamed jet emission, while in luminous 
quasar sources interactions with the ambient photon 
field will severely limit the efficiency of centrifugal 
acceleration. 

Using a simple toy model, we have shown that efli- 



cient centrifugal acceleration of electrons in the vicin- 
ity of the central supermassive black hole could nat- 
urally account for the observed VHE characteristics 
in M87, including variability on time scales of about 
two days and a hard TeV spectrum. The strength of 
the considered scenario is not only that it success- 
fully reproduces the VHE characteristics, but also 
that it allows substantial corroboration of facts by 
drawing on insight gained in accretion disk physics 
and jet formation theory. Generalizing the approach 
suggests that (i) similar to M87 other low luminous 
AGN jet sources such as, for example Gen A, may as 
well be TeV 7-ray emitting sources, and that (ii) ef- 
ficient centrifugal acceleration may lead to the onset 
of energetic particle beams in BL Lac-type objects, 
responsible for X-ray and TeV radiation on larger 
(> 10 tl) distances [27]. 
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